The transmission spectra of the graphene-based Fibonacci superlattice 
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We consider the gapped graphene superlattice (SL) constructed in accordance with the Fibonacci 
rule. Quasi-periodic modulation is due to the difference in the values of the energy gap in different 
SL elements. It is shown that the effective splitting of the allowed bands and thereby forming a 
series of gaps is realized under the normal incidence of electrons on the SL as well as under oblique 
incidence. Energy spectra reveal periodical character on the whole energy scale. The splitting of 
allowed bands is subjected to the inflation Fibonacci rule. The gap associated with the new Dirac 
point is formed in every Fibonacci generation. The location of this gap is robust against the change 
in the SL period but at the same time it is sensitive to the ratio of barrier and well widths; also it 
is weakly dependent on values of the mass term in the Hamiltonian. 
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FIG. 1 



In recent years both graphene and graphene structures 
attracted much attention which is naturally explained 
by their non-trivial properties"^. On the other hand, 
notable among the semiconductor structures including 
graphene ones are superlattices intermediate between pe- 
riodic and disordered - quasi-periodic structures, e.g. Fi- 
bonacci, Thue-Morse and others. This is due to their un- 
usual properties, such as self-similarity, fractal-like elec- 
tronic spectrum etc. Attempts has already been started 
to study some kinds of su perlat tices based on graphene 
(Fibonacci and Thue-Mors e! 10 * 11 !) which shed light on the 
behaviour of Dirac chiral fermions in the quasi-periodic 
chains. Motivated by these circumstances, in this re- 
port we study the energetical spectra of the Fibonacci SL 
based on the monolayer gapped graphene. Quasi-periodic 
modulation is due to the difference in the values of the 
mass term A in the Hamiltonian in different elements 
of the superlattice. Because of the substantial progress 
in techonology of graphene structures, in particular con- 
cerning the fabrication of the gapped graphene^^fl] we 
vary the value of A in wide range regardless of its nature 
which may be different! 1 " 2t H 

Consider the SL built of two elements "a" and "b" , see 

Fig.g 

Both elements contain the quantum well of width "w" 
and a potential barrier of width "d" . The term A a cor- 
responds to the "a" element and the barrier height is 
denoted as V a , for an clement "b" we have Af, and Vb 
respectively. The SL is constructed according to the Fi- 
bonacci inflation rule. Thus for the fifth Fibonacci gen- 



eration one can write : S f 5=abaababa. The value of the 
transmission coefficient T(E), E being the electron en- 
ergy, through the SL built for a certain generation is 
determined by the period of this generation (sequence) . 
Energy intervals in which the condition T{E) — 1 holds 
form the allowed bands in energy spectrum, gaps are as- 
sociated with values Tel. The transmission coefficient 
can be evaluated by the transfer matrix method, by ex- 
pressing T either in terms of the Green functions or the 
eigenfunctions. The latter are the solutions of the Dirac- 
like equation 
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where vf — 10 e m/s, p = (p x ,Py) the momentum oper- 
ator, a = (a x ,cr v ), a x , a y , a z — Pauli matrices for the 
pseudospin, V(x) the external potential, which depends 
only on coordinate x, I - two-dimensional unit matrix, 
the mass term is denoted by the symbol A as adopted 
in the literature. The function f is a two-component 
pseudospinor \I> = [$>a, ^b] T , ^ a, are the envelope 
functions for the graphene sublattices A and B. Sup- 
pose that the potential consists of repetitive rectangular 
barriers along the x-axis and Vj(x) = const within each 
j-th barrier. Since the y-component of the electron mo- 
mentum commutes with the Hamiltonian one can write 
*ab = ^a,b ■ e lk y' v and we get from dip : 
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where k 3 = siga(s j+ )[(E- V,-) -A 2 ] 1 / 2 , s± = E-V(x)± 
A, where units c=h = e = Vf = l are adopted. If wc 
present solutions for the eigenfunctions ^ a,b as a sum 
of plane waves traveling in the forward and backward 
directions of the x-axis then we obtain: 
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where qj = sign(s J + )y fc 2 — fc 2 if fc 2 > k 2 and qj = i ■ 
k 2 — k 2 - otherwise, = (±c/j + i ■ k y ) /k y , the top line 
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FIG. 2: Trace map for the Fibonacci SL, k y = 0, 
d = w = 0,5,V a = V b = l,A a = l,A b = 0. 




FIG. 3: Dependence of the transmission coefficient T on 
electrons energy E for the 4-th Fibonacci generation; 
values of the parameters: k y = 0, d = w = 0,5, 
V a = V b = 5, A h = 0. 



in ([3]) pertains to the sublattice A, the lower one — to the 
sublattice B. Transfer matrix connecting the wave the 
functions at thepoints x and x + Ax is found in several 
studies (see e.gP) and has the form: 




The coefficient of the transmission of quasi-particles 
through the lattice T — \t\ 2 , 



R 2 2e- l 6 ° + iW e ° - R 12 -R 2 i' w 
where 9q is the angle of incidence, and the matrix R is 

N 

expressed as a product of matrices My. R — Yl Mj, n 

i=i 

— the number of elements in the SL. 

The trace map for 6 initial Fibonacci sequences is de- 
picted in Fig. [2] for the case where we have the gapped 
graphene in elements "a" and the gapless one in elements 
"b" . First of all it is noteworthy that the quasi-periodic 
modulation by setting different values of A in different 
SL elements used in this study results in the efficient 
splitting of the allowed energy bands and thus in the 
formation of a number of gaps. And this is realized at 
normal incidence of the electrons on the lattice. 

Fig. [3] exhibits the tunneling spectrum, i.e. the depen- 
dence of the transmission coefficient T on energy E for 
the 4-th Fibonacci generation under normal incidence of 
the electron wave on the lattice; the same spectrum in 



the energy range [0, 6] is shown in Fig. [4] We see that 
the structure of some parts of the spectra is repeated 
periodically throughout the whole energy scale; one of 
these fragments of the spectrum can be considered as its 
period. The characteristic features of the period - the 
number of allowed (forbidden) bands and their widths 
change so that with E increasing the gap's widths on aver- 
age decrease. The natural result of this reduction is that 
the transmission coefficient asymptotically approaches to 
unity in the sufficiently far over-barrier region. 

The band widths don't decrease monotonically - with 
E increasing the alternate extension and reduction of 
gaps (or the allowed bands) is observed. This results 
in a wider periods (like a super-period) as it is clearly 
seen in Fig. [2] (e.g. the interval between the points L and 
P). In other words, the result of such a wavelike changes 
of band widths is the grouping of smaller structural units 
into bigger ones to form the additional structural order 
- a manifestation of the self - similarity in the problem 
considered. 

Spectra similar to that shown in Fig.|3]for the 4-th gen- 
eration are observed also for other Fibonacci sequences. 

The number of bands and the width of each of them 
depend on one hand on the SL parameters and on the 
other hand on the number of the generation. 

Here we would like to draw attention to a certain con- 
trast to the situation in conventional SL (with a parabolic 
dispersion law of the charge carriers). Unlike in the con- 
ventional SL where the calculation of bands usually is 
carried out within the barrier region, in graphene Fi- 
bonacci structures for this purpose it is advisable to 
choose the certain energy intervals, e.g. periods of spec- 
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FIG. 4: Tunneling spectrum for the 4-th Fibonacci 
generation in the energy range [0, 6]; values of the 
parameters are the same as those in Fig. [3] 



As we can see in Figs. [3j [4] at a certain energy all 
Fibonacci generations form the gap associated with the 
new Dirac point - "new Dirac point gap"^. Its location 
is almost unchanged in different Fibonacci sequences. A 
typical feature of the new Dirac point gap is that it is in- 
dep endcnt of the lattice period (d + to) and at the same 





FIG. 6 




FIG. 5 



tra (see Figs.[3j|4j) or some other fixed fragments of spec- 
tra. 

As follows from the calculations, the number of allowed 
(forbidden) bands contained in one period, in particular 
in the interval CD in Fig. [4] corresponds to the Fibonacci 
sequence and this number is subjected to the Fibonacci 
inflation rule: Z n = Z n —\ + Z n _2, where n is the number 
of the Fibonacci generation (see Fig. [2). Note that this 
rule is applied not only to the CD period but to larger pe- 
riods as well; every new super-period has its own number 
of bands. 



time is sensitive to the ratio d/w. This is proved in par- 
ticular by Fig. [5] which shows the dependence of T on E 
for the third Fibonacci generation for different values of 
w and d. Once again we draw attention to the fact that 
the band's splitting (thereby forming a series of gaps) is 
realized even in the case of the normal incidence of elec- 
trons on the lattice surface. This result is significantly 
different from what was obtained irP^ in which a quasi 

- periodic modulation was created by the difference in 
the potentials of elements "a" — the barrier and "b" - 
the well, and the splitting of bands was observed when 
k y 7^ only (oblique incidence of the wave) . 

The location of the new Dirac point depends in general 
on values of each parameter d, to, V a , Vb, A Q , A& and it 
turns out that if we put V a = Vb = V, d = w then this 
position is equal to Ed ~ ¥ and only slightly deviates 
from this value with increasing of A(V/2) is the exact 
position of the new Dirac point in the periodic lattice^) . 

The band structure for the SL of the third Fibonacci 
generation is shown in Fig. [6] in coordinates E, k y ; it 
features the spectra dependence on the angle of incidence 
9. Note that there is some extension of the new Dirac gap 
and at the same time the dependence of other gaps on 
9 is very weak (not shown in Fig. |6| . This circumstance 
is known to be common in the study of certain effects in 
graphene structures, in particular, the same result was 
stated irP: if there is a sufficiently strong effect at k y 

— then the dependence on k y is weak, see also the 
comment irP. The dashed line in Fig. [6] denotes the new 
Dirac point's location (almost equal to V/2). 
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